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Abstract 

The current distribution is studied for a finite- width two-dimensional system 
of Hall electrons, with a clear distinction drawn between the equilibrium edge 
current and the Hall current. It is pointed out that both the distribution and 
direction of the equilibrium edge current change dramatically as the number 
of electron edge states increases, and that this alternating edge current is 
another manifestation of the de Haas- van Alphen effect. The Hall-current 
distribution is substantially different from the edge current distribution, and 
it is shown numerically that the fast-traveling electrons along the sample edge 
are not the main carriers of the Hall current. 
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1. Introduction 

One basic question inherent to the quantum Hall (QHE) is where in the sample the 

current flows. The traditional explanations&^i for the QHE are based on a picture that the 
Hall current is carried by electrons in the sample bulk and regard local disorder as crucial for 
the formation of visible conductance plateaus. On the other hand, considerable attention has 
recently been directed to another picture,i^ the edge-state picture of the QHE, where the 
current is taken to flow in channels along the edges of a sample. Experiments! appear to favor 
interpretation of observed results in terms of the edge channels but little is known about 
the current distribution in real samples so far, except that some information is obtained 
from observed potential distributions.0llll Theoretically the current distribution as well as 
the potential distribution has been studied in some models, and possible connections 
between the bulk-state and edge-state pictures of the QHE have been discussed. 

The purpose of the present paper is to study the current-carrying properties of the 
electrons in the sample bulk and near the edges in detail. In equilibrium a prominent current 
flow arises along the edges of a Hall sample. It is pointed out that both the distribution 
and direction of the edge current change dramatically as electron edge states increase in 
number and that this alternating edge current is another manifestation of the de Haas - 
van Alphen (dHvA) effect. The Hall-current distribution turns out substantially different 
from the equilibrium current distribution, and it is shown numerically that the fast-traveling 
electrons along the sample edges, though carrying a large amount of current per electron, 
are not the main carriers of the Hall current. 

In Sec. 2 we study the equilibrium current distribution in a disorder-free Hall sample of 
flnite width. In Sec. 3 we examine the Hall current distribution and discuss the transport 
properties of the edge current. Sec. 4 is devoted to concluding remarks where the effect of 
disorder is discussed. 
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2. Current distributions in equilibrium 

Consider electrons confined to an infinitely long strip of width W (or formally, a strip 
bent into a loop of circumference ^ W) in the presence of a uniform magnetic field B 
normal to the plane, described by the Hamiltonian: 

H, = ^co{epl + il/i')iy-yor}, (2.1) 

written in terms oi lu = eB/m, the magnetic length i = 1/V eB and yo = px/{eB); the 
Landau-gauge vector potential {—By, 0) has been used to supply a uniform magnetic field. 
We take explicit account of the two edges y = ±W/2, where the wave function is bound to 
vanish. 

The eigenstates of Hq in the sample bulk are Landau levels with discrete spectrum 
Lj{n + |) labeled by integers n = 0, 1, 2, ■ ■ -, and degenerate in yQ = i'^Px', the yo labels the 
center-of-mass position of each orbiting electron. The eigenfunctions in the presence of sharp 
edges are still labeled by = (n, yo): 

V'^(a;,y) = (L,.)-^/'e"^^0^(y), (2.2) 

with (pNiy) giveni by the parabolic cylinder functions^ Dy{-^\/2{y — yo)/^) for electrons 
residing near the edges y = The energy spectrum of each Landau level n 

enivo) = uj |z/n(yo) + ^| (2.3) 

is fixed from the boundary condition (pN = at y = ±W/2. The spectra Vniyo) are deter- 
mined numerically and are functions of a dimensionless combination y^ jt- 

yn{y,) = yn[y^/^i (2.4) 

where y^ = yo + W/2 refers to the value of yo measured relative to yo = —W/2 near the 
"lower" edge; = yo — W/2 near the "upper" edge yo ~ W/2. As shown in Fig. 1, 
^n{yo) = ^n{—yo) rise sharply only in the vicinity of the edges yo ~ ±W/2, and recover the 
bulk values n as yo moves to the interior a few magnetic lengths away from the edges;il for 
the lowest n = level, < z/g < 0.003 for y^ > 2.5i. Applying the Feynman-Hellman 
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theorem to dHo/dyo reveals that a single electron state N — {n,yo) has its center-of-mass 
at position 

yl-^^{N\y\N)=yo-fun'{yo) (2.5) 

on the y axis. This y^ stays within the interval [—W/2,W/2] while yo — Px can take 
values outside it for electrons residing near the edges. 

The normalized wave functions (pNiy), taken to be real, are functions of y — yo and y^: 

My) = (i/v^)0n ((y - yo) A; y^/i) , (2.6) 

where we have taken 0„ dimensionless so that its dependence on dimensional quantities i, yo, 
etc. is made exphcit. Let us denote the spatial distribution of a single electron state N as 

1 2 

PN(y) = pn{y; yo) = j^n {{y - yo)/^; y7 1^) > (2-7) 

which is uniform in x and which is normalized so that 

/ dypniy;yo)^l. (2.8) 

J-W/2 

The associated current density per electron is written as a product of charge (— e), velocity 

{p^ - eBy)/m and density Pn/ L^: 

My\N) = {y-y^)p^{y;y^). (2.9) 
^x 

The distributions p„(y; yo) are highly localized and are nonzero only for \y — yol < 0{i). 
In the sample bulk [i.e., |yo| < W/2 — 0{£)], in particular, p„ are even functions of y — yo: 

pT\y - yo) = H^{Y)' e-^V(£2'^n!0F), (2.10) 

with y = (y — yo)/£. The profile of Pn{y]yo) gets gradually deformed as yo lies closer to 
the edge, as shown in Fig. 2 (a) for the n = level. There the y^ — 3.5£ electron state, 
depicted with a solid curve, has uq ^ 10^^ and thus exhibits essential characteristics of an 
electron in the bulk; clearly the associated current distribution, depicted with a solid curve 
in Fig. 2 (b), is characteristic of an orbiting electron and carries no net current. In contrast, 
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for an electron lying closer to the edge the cyclotron motion gets deformed, yielding a net 



The current-carrying property of each electron state is made manifest by decomposing 
the current jx{y\N) into the circulating component oc {y — y'^) Pn{y',yo) (that carries no 
net current) and the traveling component oc —i'^h'n{yo)Pn{y',yo)', see Fig. 2 (c) and (d). It 
is clearly seen from these figures that the electrons residing near the edge (the edge states) 
carry a large amount of current. Classically these edge states are visualized as electrons 
hopping along the sample edges.i^ They travel fast (with velocity = ujPun ~ ±uji) in 
opposite directions at opposite sample edges. 

Let us next examine the current distributions filled Landau levels support. We shall 
focus on the case of well-filled levels close to integer filling. It is convenient to use the Fermi 
potential Ep = uj{np + 1/2) or rip to specify the filling of each level: For a given rip Landau 
levels n with n < np are populated densely with electron states over the range (i/(7)n < Ho < 
(yo)n so that t'n((z/o^)n) = np] the filling fraction of each level is /„ = [iyo)n~{yo)n]/W with 
the number of electrons given by {Ne)n = i^B / 2n) L^W fn- The np stays constant np = n 
while the nth level is gradually filled with the electron states residing in the sample bulk 
(the bulk states), and it increases as soon as the edge states start to be filled. Such filled 
states (n, yo) give rise to the current distribution 



In view of Eq. ( p.lOp , it is clear that jx{y) = in the sample bulk. Figure 3(a) shows the 
current distribution near the lower edge y = —W/2 for various values of np. For < np < 1 
only the n = level carries current, for 1 < np < 2 two levels [n = 0, 1) contribute, and so 
on. 

There are two notable features read from the figure: (1) The current flows along the 
edge in a channel whose width ~ 0{i) increases with the number of current-carrying levels. 
This feature has been well knownSi and is naturally expected from the current fraction 



current which, in view of Eq. ( |2.5| ), is proportional to the gradient of the spectra i/„: 




(2.11) 




(2.12) 
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per electron in Eq. ( p.ll|) . Somewhat unexpected is the following: (2) Both the spatial 



distribution and direction of the current change dramatically with rip. For small rip ~ 0.1 
the edge current flows predominantly in the negative x direction while a large fraction of 
current flows in the opposite direction for rip ~ 0.9; the integrated amount of current changes 
its sign at np = 0.5. Analogous patterns of current flow are also seen for each integer interval 
n < rip < n + 1. The amount of current integrated over one edge oscillates with np, as shown 
in Fig. 3(b). 

These changing and oscillating current distributions may appear puzzling if one notes 
Eq. alone. There is, however, a simple resolutionil: Associated with each orbiting 



electron is a circulating current oc {y — y'^) PnijJ] yo) that produces magnetic moment point- 
ing opposite to the applied magnetic field B. For a dense collection of electrons in the sample 
interior this circulating diamagnetic current is averaged out to vanish locally, but leaves a 
current circulating along its periphery. Indeed, collecting the circulating current components 
alone yields a current distribution j™'^{y) localized near the edge, as shown in Fig. 4 (a) for 
the n = level with different np. Also shown in Fig. 4 (b) is the current distribution j^^{y) 
formed by the traveling current components alone for the same n = level. They combine to 
build up the current distributions of Fig. 3(a). It is now clear that the current distribution 
for = 0.1 in Fig. 3 (a) primarily derives from the orbital diamagnetic current whereas the 
edge-driven traveling current dominates in the np = 0.9 case. Note that the edge-driven cur- 
rent works to cancel the orbital diamagnetic current. The cancellation is partial in quantum 
theory,0 leading to the Landau diamagnetism of electrons for three-dimensional samples. 

In the present two-dimensional case it is possible to calculate both the diamagnetic and 
edge-driven components of the edge current explicitly. As a preliminary step, let us first 
consider a collection of electron bulk states (of the nth level) that fill up a half-infinite 
interval yo > 0. They lead to the current distribution 



jr'(z/) = y;^' ^^^pr (^) (2.13) 

localized around y = with spread Ay = 0{i); for n = this is given by — (2y^)"^e"^^/^^ 
times {eu/2TTi). The integrated amount of this diamagnetic current increases with n: 



/ dyjr\y) = -(n + l)'^. (2.14) 

Consider next a collection of electron states that fill up the interval {yo)n ^ Vo ^ V with rj 
lying somewhere far in the sample interior. The circulating current components associated 
with these states, though carrying no net current, combine to build up two prominent 
current distributions localized near the edge y = —W/2 and around y = rj, which are 
equal in net amount of current and opposite in sign. This shows that the diamagnetic 
component of the current localized near the edge y = —W/2 supports a fixed amountH 
equal to — (ec<j/27r)(n + 1/2), irrespective of its distribution as well as the shape of edge 
potentials. In particular, the four different distributions of j'^^'^iy) in Fig. 4(a) support the 
same amount of current equal to — (eti;/47r). 

On the other hand, the edge states of the rath level arise for n-p > n and, as seen 
from Eq. ( |2.11| ), carry the amount of current equal to {euo /2Ti){n-p — n) at the lower edge 



y ~ —W/2. With the two current components put together the nth level alone supports the 
amount of edge current 

= — (np - 2n - 1/2) fornp > n. (2.15) 

27r 

For the integer interval n < n-p < n + 1 the lower (n + 1) filled levels combine to support 
the amount of current 



J-- = ^(n + l) 



np-{n + 



(2.16) 



localized near the edge y = —W/2. This current-np relation accounts for the alternating 
edge current of Fig. 3(b), and shows that the (integrated) edge current changes its direction 
precisely at np = n + 1/2. 

The alternating edge current of Fig. 3 is intimately connected to an oscillatory variation of 
magnetic susceptibility of an electron gas with varying magnetic field, known as the de Haas 
- van Alphen (dHvA) effect.i^'ll Indeed, the current J'^'^s^^ circulating along the sample edge, 
gives rise to uniform magnetization normal to the sample plane and of magnitude (per unit 



area 
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= ^B{l/27if){n + l)[2(nF - ^) - 1] (2.17) 

for n < Up < n + 1, apart from corrections of 0{i/W) that vanish as W ^ oo, where 
Hb = euj£'^/2 = e/{2m) is the Bohr magneton (with effective mass m). 

This result is readily verified by thermodynamic methods as well. The simplest way is 
to consider how the energy of the present finite-width system, 

E = iLj27Tf) Y^fdyoiu {i^nK/i] + 1/2} , (2.18) 

responds to an infinitesimal variation of the magnetic field B. Keep = VqI^ fixed and 
calculate the magnetization = —{l/WLx)SE/6B: 

= / ^ {2^n(yo) + 1 - {2yo - y^yniVo)} . (2.19) 

Note that i/„ = n in the sample bulk and that i^n ^ only in the edge regions. The yo integral 
therefore is equal to [2n + 1 — 2 (rip — n)] for each filled level n, apart from corrections of 
0{i/W); this result is consistent with Eq. ( |2.15| ) and thus recovers Eq. ( p.l7| ). It is also 
enlightening to confirm Eq. ( p.l9[ ) by a direct calculation of the thermodynamic potential 
T; B) = -T E» ln(l + e^'^-^')/'^^) with T ^ 0. 

With the spectra z/„(?/o) determined numerically, it is a simple task to express the mag- 
netization ( p.l9| ) as a function of B or the filling factor / = J2n fn oc N^/B. As seen from 
Fig. 5, the magnetization per electron, = Mz/{Ne/WLx), plotted as a function of / 
shows an oscillatory variation characteristic of the dHvA effect. As usual, the gradual varia- 
tion of niz in each integer interval n < f < n+1 — 0{i/W) is ascribed to the diamagnetism 
of orbiting electrons in the bulk. In each tiny interval n + 1 — 0{i/W) < / < n + 1, 
makes a drastic change from diamagnetism to paramagnetism that is caused by the edge 
states of the nth level. 

Note here that the electron edge states, because of their tiny filling fraction~ 0{i/W), 
scarcely contribute to the internal energy E. Thus, even without explicit account of the edge 
states, one can still arrive at the dHvA effect by thermodynamic methods (i.e., through E or 
fl), except that now shows discontinuous jumps at integer filling f = n. The edge states 
are invisible in such derivations but are certainly present physically: Note that magnetization 
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by nature is continuous as a function of B/Ne-, as seen clearly from the W = 30i case of 
Fig. 5. In view of this continuity, the prominent sign change of near integer filling f = n 
(which can be derived by thermodynamics without the edge states) does imply the presence 
of the edge states and, in particular, the alternation in direction of the edge current. In this 
sense, the alternating edge current in Fig. 3 is another aspect of the dHvA effect. 

The alternation of the edge current is caused by competition between the circulating 
and edge-driven components of the current near the edge. The edge-driven component is 
naturally very sensitive to the shape of the edge potential while the circulating component 
is not; see Fig. 2 (c) and (d). Correspondingly the edge current in general varies in distribu- 
tion according to the shape of edge potentials. Still, in net amount per edge, both current 
components remain unaltered so that Eqs. (|2.15|) and (|2.16|) hold, irrespective of the details 
of edge potentials. 

3. Hall-current distributions 

The alternation of the edge current, which we have just seen, does not necessarily imply 
that the information the edge current carries also changes the direction of propagation. In 
this section we clarify this point by studying how the edge current responds to a Hall field. 

For simplicity let us consider a uniform field Ey = —dyAo{y). Its effect is readily taken 
care of by making the shift 



in Hq of Eq. ( |2.1| ), and the full Hamiltonian Hq — eAoly) leads to the new spectrum 
^n{yo) + eEyyo + 0{Ey). The normalized wave function 4>N{y) is given by Eq. ( p.6|) 
with the replacement yo yo, and the current per electron jx{y\N) by Eq. ( |2.9| ) with 
Pn{y',yo) Pn{y',yo)- The current distribution jx{y\N) itself barely changes thereby be- 
cause, under conditions of practical interest, the deviation yo — yo is negligibly small. 
Actually the Hall current we are interested in is the small deviation. 



yo^yo = yo- {ef/uj)E, 



'y 



(3.1) 



Nx{y\N) = — {y- yo) {pn{y; yo) - pn{y; yo)} 



(3.2) 
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representing the response to an applied field. As seen from Fig. 6(a), unlike jx{y\N), the Hall 
current per electron Ajj.{y\N) is primarily composed of traveling components. Here we see 
that, while the edge states and bulk states are drastically different in the amount of current 
they carry, they are essentially the same in the Hall-current transport. Numerically an edge 
state supports even a smaller amount of Hall current than a bulk state, as seen from Fig. 6(b), 
where the numerically- integrated amount of Hall current per electron, AJ^ = J dy Ajx{y\N), 
is compared with the net current per electron, Jx = J dy jx{y\N) oc —Uniyo)- 

It is possible to understand the current-carrying properties of each electron state in a 
more general way: Of the current jx{y\N), the circulating component oc (y — y^™) p„ carries 
no net current. It is the traveling component, associated with the drift of a Hall electron 
with velocity 

c^^' Vn{yo) + Ey/B, (3.3) 

that carries a net current (and hence information with it). Accordingly disturbances caused 
upon a sample, e.g., by varying a magnetic field or electron population will propagate in 
a direction fixed by the edge with velocity~ ujI. In contrast, the effect of a Hall field 
(i.e., the Hall current) propagates in a direction fixed by the polarity of Ey with velocity 
~ Ey/B. [Numerically ~ 10^ cm/s for typical values ~ 10 meV and I ~ 100 A while 
Ey/B ~ 10^ cm/s for Ey = I V/cm and 5 = 5 T.] It follows from Eq. (jSj) that the net 
amount of Hall current per electron is proportional to 1 — £^ Wiyo), which reproduces the 
numerical result in Fig. 6(b) very accurately. 

Filled Landau levels support the prominent edge-current distributions of Fig. 3(a), which 
remain essentially unchanged in the presence of a Hall field as well. These edge currents flow 
in opposite directions at the two opposite edges of a sample, and in equilibrium with Ey = 
they combine to vanish.ili When a Hall field Ey is turned on, the Hall current emerges as 
a small deviation from the equilibrium distribution, as shown in Fig. 6(c) for the present 
impurity-free and uniform- E'j^ case. It is clear from Fig. 6, contrary to some tempting ideal, 
that the electron edge states, because of their tiny filling fraction ~ 0{i/W), share only a 
tiny portion of the Hall current .0 
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4. Concluding remarks 

In this paper we liave examined current distributions in a Hall bar in the regime of 
the integer QHE, and shown, in particular, that the edge current changes its distribution 
and direction as the number of electron edge states increases. This dramatic change is 
a consequence of competition between the circulating diamagnetic component and edge- 
driven traveling component of the current carried by electrons near the sample edge, and 
is closely related to the dHvA effect. It should be emphasized that the dHvA oscillation 
of magnetization is indirect evidence for the edge states; it is clear physically that without 
them no oscillation would arise, since orbital magnetization alone leads to diamagnetism.0 
In this connection we have seen explicitly that magnetization is continuous as a function of 
B /Ne when the edge states are properly taken into account. We have also seen that the Hall 
current flows in a manner quite independent of the equilibrium edge-current distribution. 

Finally we would like to comment on the influence of disorder. Rapid motions of Hall elec- 
trons like cyclotron motion and acceleration by the edge are potentially not very sensitive to 
disorder. The edge current of Fig. 3, resulting from such rapid motions, will therefore remain 
prominent in the presence of weak disorder as well, and continue to alternate in direction 
as the number of edge states increases. In contrast, disorder will substantially modify the 
current carried by slowly-drifting electrons in the sample bulk: In the presence of disorder 
a large fraction of electron bulk states get localized and cease to carry current; at the same 
time, the surviving extended states carry more current and achieve exact compensation.!'! 
This exact current compensation is a consequence of electromagnetic gauge invariance and 
takes place under general circumstances involving both bulk and edge states.0 

Note now that, of the electron bulk states, those residing near the edge of the sample 
bulk ("bulk edge"), though influenced by disorder, would have a better chance of staying 
extended than those far in the bulk. In view of current compensation, a considerable portion 
of the Hall current, redistributed via disorder, would therefore flow along the sample edges. 
In this way it has been pointed out0 that there are two kinds, fast and slow, of edge current 
in Hall samples in the regime of the integer QHE. 
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These two kinds of edge current differ in channel width and in direction of flow. The 
fast component is nothing but the alternating edge current discussed so far, consisting of 
the two (circulating and traveling) subcomponents; it has a channel width of 0{i) and 
flows in opposite directions at opposite sample edges. In contrast, the slow component, 
the "bulk-edge" Hall current, will have a channel width, related to some localization length 
characteristic to the bulk edge, which could well be larger than 0{i). This slow edge current 
will flow in the same direction at opposite edges and reverse direction when the polarity of 
the Hall field is flipped; an observed Hall-potential distribution0 appears to be in support 
of this feature. 

A numerical experiment is now under way to study the current distribution for small 
samples with randomly distributed impurities, and is yielding results that appear to confirm 
the effect of disorder on the current distribution described above; details will be reported 
elsewhere. 
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FIGURES 
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FIG. 2. (a) Spatial distributions PniViVo) of single electron states with different i^o(yo) in 
the n = level. The solid curve represents the « 10~^ electron state with y^^ = 3.5^, the 
long-dashed curve the vq = 0.1 state with 1.40^, and the dashed curve the vq = 0.5 state 

with 2/0^ 0.541^. The coordinate = y + W/2, defined relative to the edge y = —W/2, 
is measured in units of i. (b) Distributions of the associated current per electron. Here, for 
convenience, the current jxi,y\-^^ is plotted in units of — euj/Lx so that its sign refers to that of 
the velocity Vx along the x axis, (c) Decomposition of jx{y\N) into the circulating component 
jx^^{y\N) and the traveling component j^{y\N). 
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FIG. 3. (a) Current distributions filled Landau levels support near the sample edge y = —W/2; 

np = Ey/lo — 1/2 refers to the Fermi potential Ep. The current jxiy) is measured in units of 
—eoj/{2iT) and the coordinate = y + W/2 in units of £. (b) The integrated amount of the 
current per edge [in units of — ea;/(27r)] oscillates with np. 
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FIG. 4. (a) Current distributions near the edge y = —W/2 built up by the circulating current 
components alone for the n = level with different values of np. The np ~ 10"^ case refers to the 
n = level filled with electron states > 3.M. The current is measured in units of — ea;/(27r), 
and the coordinate y^ in units of £. (b) Current distributions formed by the traveling current 
components alone for the same n = level. 
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FIG. 5. Magnetization per electron [in units of v.s. the filling factor / for a Hall bar 

of finite width W = 30^, 50^, and 100^. 
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FIG. 6. (a) Distributions of the Hall current per electron, Ajx{y\N), for electron states 
with different fo(yo) in the n = level. The current is plotted in units of 10~^ x {—eu/Lx) 
for e£Ey/uj = 10~^. (b) The Hall current per electron AJ^ = / dy Ajx{y\N) (in units of 
—e^i'^Ey/Lx) for electron states in the n = level, in comparoson with the net current per electron 
Jx = J dyjx{y\N) (in units of euj£/Lx). (c) Hall-current distributions near the edge for the n = 
level. The edge states increase in number with np. 
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